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I. INTRODUCTION 



At low energies, quantum chromodynamics can be described by a low-energy effective 
theory, chiral perturbation theory %PT. To leading order, it has two a priori unknown 
parameters, the chiral condensate and the pion decay constant. In this paper, we compute 
these constants by an ab initio computation. Our strategy itself is not new: we compute 
two-point functions in the epsilon regime and fit them to their %PT predictions. However, 
for chiral perturbation theory to be an accurate description of low-energy phenomena, we 
need calculations on the QCD side at very low pion masses, probably lower than typically 
reached in current simulations, and at the same time on large enough volumes to control 
higher order corrections. 

Basically, there are two regions in parameter space where calculation of chiral pertur- 
bation theory can be carried out: one is the so-called p-regime where we are essentially at 
infinite volume. The pion wave length is much smaller than the size of the box, and the 
small corrections due to its finite extent can be taken into account analytically. 

The other region is the e-regimeQ. There the pion wave length is much larger than the 
spatial and temporal size of the lattice. One can carry out the integrals over the (effectively) 
constant pion modes exactly and arrive at a different power counting. This regime has 
considerable appeal. First of all, even at next-to-leading order, only the two leading order 
low-energy constants of chiral perturbation theory enter into observablesjjjl, Q|. Second, 
being fundamentally a finite volume regime, once one has reached it, lowering the pion mass 
actually improves the validity of x?T at a given order, whereas in the p-regime one needs 
to increase the volume at the same time. 

The problem, however, lies in the task to actually reach this regime in numerical simula- 
tions. Because the epsilon expansion is in powers of 1/(FL) 2 , it turns out that we actually 
need a fairly big volume and therefore a very small quark mass to satisfy the condition 
m n < 1/L. This poses many algorithmical problems. In a recent publication \^ we have 
presented a setup with which we can actually reach the regime of small quark masses at 
moderate cost avoiding many of the problems a more direct approach would face. We pro- 
pose to use smeared link improved Wilson fermions 0, 0] to generate an ensemble of gauge 
configurations above the desired quark mass. From this ensemble, we reweight to the quark 
mass which we actually want to reach. Using this method, we manage to reduce the sea 
quark mass by roughly a factor of 2 to 4. 

Reweighting has another advantage apart from the actual possibility to go to such light 
quarks. At very small quark mass statistical fluctuations grow dramatically because the 
mass no longer provides an infrared cutoff. Therefore, very small eigenvalues of the Dirac 
operator can appear, which are suppressed by the fermion determinant, but also lead to large 
values of the measured observables. We thus have an anticorrelation between the weight and 
the function value which means importance sampling breaks down. Reweighting avoids this 
problem. The small eigenvalues are not as efficiently suppressed, the region of large signal 
gets over-sampled and the error is actually reduced. 

One might wonder how the explicit chiral symmetry breaking of Wilson fermions influence 
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simulations done with very light quarks in the epsilon regime. Unlike in the p-regime, the 
finite volume of the e-regime provides the system with an effective, dynamical IR cut-off. 
The chiral symmetry violations will remain managable even in the chiral limit as long as 
they are small compared to the inverse lattice size. The numerical study of the continuum 
limit to prove this expectation will be the subject of a future work. 

In the current paper, we apply reweighting to a large volume data set and provide mea- 
surements of the low-energy constants F and S. We already mentioned that this t ype of sim- 
ulation is not new. A number of quenched studies have been carried out 
using overlap fermions, which proved the feasibility of the method but also highlighted the 
need for sufficie ntly large volume. Recently, there also have been computations with dy- 
namical overlap jljj and twisted mass 15fl fermions to which we will compare our results in 



the conclusion. For a recent review of these calculations see Ref. [16|. The strength of our 
calculation is that we compare two different volumes, of which the larger one has not been 
reached in previous studies. The simulations in this work were fairly inexpensive, and it is 
within reach to repeat the calculation at a finer lattice spacing to verify scaling. 

The outline of this paper is the following: We first describe in Sec. [Til the set-up of our 
simulations, the generation of the ensemble of gauge configurations and the details of our 
reweighting procedure. The relevant renormalization constants are computed in Sec. DTD 
using the RI-MOM scheme. In Sec. [TV] we collect the e-regime formulas, the procedure of 
the extraction of the low-energy constants and give the results. 



II. SIMULATIONS 

The numerical simulations for this project were done with 2 flavors of nHYP smeared 
Wilson-clover fermions and one-loop Symanzik improved gauge action. The action and 
the simulation method are described in details in Refs. 0, 0|. We use tree-level csw = 
1.0 clover coefficient, so our action is not fully O(a) improved. Based on our quenched 
investigation we expect that a nonperturbatively improved action would require csw ^ 
1.2, so even with only tree-level improvement the 0(a) corrections are likely small. We have 
generated two sets of gauge ensembles, both at gauge coupling (3 = 7.2. The first set consists 
of 16 4 configurations at n = 0.1278, the second 24 4 configurations at n = 0.12805. We have 
180 and 154 thermalized configurations, separated by 5 trajectories at the two volumes. The 
autocorrelation is 3-4 trajectories for the plaquette, and about the same for the two point 
functions. Preliminary results using the first set were already reported in Ref. 

We set the lattice scale from the static quark potential, using r = 0.49 fm for the Sommer 
parameter. On both configuration sets we found r /a = 4.25(2) (the error is from the larger 
volume set where we have a better signal for the potential), giving a = 0.1153(5) fm. With 
this value the physical volumes are (1.85 fm) 4 and (2.77 fm) 4 . Based on the PCAC quark 
mass and the pseudoscalar and axialvector renormalization factors (see Sec. lIIII) . we estimate 
the renormalized quark mass in the MS scheme at 2 GeV to be 22 and 8.5MeV, respectively. 
These values, and some other details of the simulation, are listed in Table fl] 

The dynamical simulations were performed at particularly low quark masses, even if we 
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180 
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0.1279 
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180 


0.0088(5) 


16.5 




0.1280 
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180 


0.0058(7) 
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0.12805 


16 


180 


0.0047(8) 
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0.1281 


16 


180 


0.0028(11) 
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0.12805 


0.12805 


24 


154 


0.0044(3) 


8.5 




0.12810 


24 


154 


0.0030(3) 


5.8 




0.128125 


24 


154 


0.0024(3) 


4.2 




0.12815 


24 


154 


0.0019(4) 


3.8 



Table I: The parameters of the simulation. The first column gives the coupling k of the dynamical 
simulation, the second the reweighted coupling K Tew . The last column is the renormalized quark 
mass using m = mpcAC^/^p- 



consider the relatively large volumes. This is possible due to the highly improved chiral 
properties of the nHYP smeared clover fermions. Figure Q] shows the histogram of the 
absolute value of the lowest Hermitian eigenmode of the configurations. Simulations with 
Wilson-like fermions require a well defined gap between zero and the first eigenmode 
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As the figure shows, our simulations are safe on both volumes, though very close to the low 
mass limit. In both cases, the median of the distribution is about four times its width o", 
which is above the 3cr stability criterion of Ref. [18|. However, it also shows that going lower 
in the quark mass can be very dangerous because of the algorithm becoming unstable. The 
same paper predicts VVa/a to be a scaling quantity. We measure 0.56(4) and 0.77(5) for 
the 24 4 and 16 4 ensembles respectively. 

Starting from the original configurations one can explore a range of quark masses in fully 
dynamical systems by reweighting the configurations. In Ref. Q we have described an effec- 
tive technique to calculate the necessary weight factors. It is a stochastic calculation, and 
one must take care not to introduce significant statistical errors with the stochastic process. 
We apply three methods, low mode separation, determinant breakup, and ultraviolet (UV) 
noise reduction to control the statistical fluctuations. In both the 16 4 and 24 4 ensembles 
we separate 6 low Hermitian eigenmodes. In addition we break up the determinant to the 
product of 33 and 60 terms for each Ak = 0.0001 shift in reweighting on the 16 4 and 24 4 
volumes, respectively. To control and remove some of the UV noise we introduce a pure 
gauge term in the reweighted action . This term is just an nHYP plaquette term and has a 
very small coefficient. We found that in our system it can be chosen to be proportional to 
the shift in k, 

AiHYP = 6.0(/t - K rcw ) (1) 

on both the 16 4 and 24 4 configuration sets. This value is so small that there is no difference 
within the errors in the lattice spacings or quark masses between /3 n HYP = and Equation 
[U For further details and the exact definition of the reweighting action we refer to Ref. 
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Figure 1: The Hermitian gap distribution on the original 16 and 24 4 configurations. The dashed 
lines correspond to ampcAc • 
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Figure 2: The distribution of the pseudoscalar correlator at t = 8 on the 16 4 ensemble at k = 0.1280. 
On both panels the dashed line is the partially quenched distribution and the solid lines correspond 
to the reweighted distributions. Left panel: /3 n HYP = 0, right panel: /3 n HYP as in Equation [TJ 

especially Equation 17 and Figure 4. 

In addition to removing the UV fluctuations, the introduction of the nHYP plaquette 
term also increases the overlap between the original and target ensembles. The largest 
weight factors are pushed from the edge of the plaquette distribution to the middle, where 
the statistical sampling is better, and the effect is similar for other observables as well. We 
illustrate this in Figure [2] with the distribution of the pseudoscalar correlator at t — 8. All 
data correspond to the 16 4 data set at n = 0.1280. The left panel shows the reweighted 
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Figure 3: The Hermitian gap distribution on the original and lightest reweighted ensembles for 
both volumes. The histograms are labeled by the corresponding lattice quark mass which is also 
indicated by the dashed lines. 

distribution without the nHYP plaquette term, the right one with /3 n HYP as in Equation CD 
For reference both panels show the partially quenched (unweighted) distribution. In both 
cases the overlap between the reweighted and partially quenched distributions is excellent, 
there is no sign that reweighting would prefer region that is poorly sampled by the original 
ensemble. The main difference between the reweighted and partially quenched data is the 
suppression of the long tail of the latter one, a quenching artifact. The apparent increase 
of the reweighted distributions is mainly due to normalization: the dynamical distribution 
is narrower, resulting in a higher peak. It is worthwhile to emphasize that including the 
nHYP plaquette term does not introduce any systematic error, rather it improves the overlap 
between the ensembles, especially at larger mass difference. 

Even though there is no strong difference between the two panels of Figure the in- 
troduction of the nHYP plaquette term reduces the statistical errors by up to 40% for our 
lightest 16 4 data set, and the results we present in the following reflect that. On the other 
hand, on the 24 4 volumes within the k range we reweight to there is no difference between 
the actions with or without the nHYP plaquette term, and the results we present here were 
obtained with /? n HYP = . 

With reweighting it is possible to reach an eigenvalue that is negative or at least smaller 
than the typical eigenvalues of the Dirac operator. We approach that case with our last 
coupling on the 16 4 ensemble where at least one configurations has a negative real Dirac 
eigenvalue and several has nearly zero eigenvalues, and even more on the lightest reweighted 
24 4 configurations where 5% of the configurations have a negative Dirac eigenvalue and 
even more have a nearly zero one. These configurations are suppressed by the weight factor, 
nevertheless as we will see later one encounters increased statistical errors on these ensembles. 
In Figure[3]we compare the Hermitian gap distribution on the original and lightest reweighted 
ensembles for both volumes. The distribution shifts towards zero but configurations with 
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negative or near-zero eigenmodes are strongly suppressed, and that maintains the gap. While 
the PCAC quark mass approches zero, the median of the gap, controlled by the finite volume, 
remains finite. 



III. RENORMALIZATION FACTORS 

In order to connect the lattice meson correlators to the physical ones we have to deter- 
mine the corresponding renormalization factors. We used the standard RI-MOM method 



191 ]. where one calculates bilinear quark operators (p\Or\p) at specific lattice momentum 



p = n and matches them to the corresponding tree-level matrix element. Afterwards 



20, 21 . The 



the lattice values are connected to the continuum MS scheme perturbatively 
renormalization scale \i has to be much smaller than the lattice cut-off to minimize lattice 
artifacts but much larger than the QCD scale for continuum perturbation theory to work 
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Our code to calculate the lattice matrix elements is based on the one used in Ref 
We used 80 propagators from the 16 4 data set to calculate the vector, axialvector, scalar and 
pseudoscalar matching factors in the chiral limit. While most dynamical calculations do the 



chiral extrapolations on partial quenched data [23|],[2J], we can do this extrapolation on fully 



dynamical configurations on our reweighted ensembles. We used 5 k values, 0.1278, 0.12785, 
0.1279, 0.12795 and 0.1280, corresponding to quark masses 10-20 MeV. We extrapolated the 
vector, axial and scalar data linearly in the quark mass, though the data shows no mass 
dependence within errors. This is not that surprising, since our quark masses are light. The 
pseudoscalar density couples to the Goldstone boson channel and it develops an 0{l/m) 
singularity in the chiral limit. We subtract this pole assuming a linear mass dependence 
for the quantity m/Zp. Again, with our light mass values we expect this assumption to 



25, 26 . Nevertheless the 



hold, and our data are indeed consistent with a linear dependence 
subtraction introduces fairly large errors at small /i. 

In Figure [4] we show all four renormalization factors converted to the MS scheme at 
2GeV, as the function of the original lattice momentum p 2 = fi 2 . The vector and axialvector 
factors are scale independent, any deviation from a constant is due to lattice artifacts. In 
our case Z v is constant over the whole range, while Za shows a slight drift at larger /i values. 
Calculations with Wilson and Wilson-like improved fermions show similar trends for these 



quantities [23, ]25|, [26 



The scalar and pseudoscalar operators depend on the energy scale. We connect the lattice 
data to the continuum one at identical energy, then, using the known 3-loop expression for 
the running of the coupling, run it to fi = 2GeV . We plot these values, therefore Zp and 
Z s in Figure [4] should also be constant. Because of the subtraction of the Goldstone pole 
the errors are large at small fi for the pseudoscalar, but we find a long, stable plateau at 
larger lattice scale values. The scalar operator, on the other hand, shows quite large lattice 
artifacts at higher scales. Again, this trend has been observed before with other actions. 
The horizontal lines in Figure [4] indicate the range where we extract the renormalization 
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Figure 4: The renormalization factors for the vector, axialvector, pseudoscalar and scalar operators 
as the function of the lattice energy scale. All values are converted to the continuum MS scheme 
at \x = 2GeV. 



factors. Our final values are 



yMS 



Z^ s (2GeV) 
Zlp(2GeV) 



0.96(1) 
0.99(2) 
0.90(2) 
1.01(3) 



(2) 



As a simple check we compare the renormalized quark mass as predicted form the bare quark 
mass m = Z^ x m hl 1/(2k) — 1/ (2/€ cr ) , and from the PCAC mass m r = Z^Z^mpcAc- Fitting 
^pcac linearly in 1/(2k) we predict n cr = 0.12821 and from the slope ZpZ^Z^ 1 = 0.94(3). 
This is consistent from the value obtained from Equation [2, 0.90(7). The fact that all four 
matching factors are close to one indicates small perturbative corrections, as it is usually 
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seen with smeared link actions. 



IV. e- REGIME ANALYSIS 

In the e-regime the pion correlation length is large compared to the linear size of the 
lattice, the light pseudoscalar mesons dominate the dynamics. Nevertheless, in order to 
incorporate the massive modes the volume has to be large compared to the QCD scale. One 
assumes that the quark mass is m = C(e 4 ) and the inverse size 1/L — 0(e) ( L 4 = V = 
LgL t ), but 1/L 3> Aqcd- The dimensionless quantity mEV, or equivalently m^F^V , is kept 
order one. Chiral perturbation theory predictions are organized in power of e 2 or 1/(FL) 2 . 
Predictions for various meson correlators are known to next-to-leading (NLO) order (C(e 4 )), 
except for the pseudoscalar that has been calculated up to C(e 6 ). In our fits we use the 
NLO predictions as those depends only on two low-energy constants, £ = \im m ^ Q (qq) and 

F = lim m ^ F w a ana. 

These xPT results are based on a chiral (continuum) action. One 
expects extra terms, due to the explicit chiral symmetry violation of the Wilson fermion 
action, in our situation. However these corrections typically show up at the same order as 
the higher order chiral constants L 3 ,L^, i.e. only at next-to-next-to-leading order in the 
epsilon regime. 

As the quark mass decreases in a large volume (p-regime) simulation, the chiral symmetry 
breaking effects of Wilson fermions get large compared to the mass, and that can create large 
lattice artifacts. In practice the continuum limit has to be taken before the chiral limit. The 
situation is different in the e-regime, where the finite volume of the system creates an infrared 
cutoff even at vanishing quark mass. This effect is well illustrated by the Hermitian gap 
distribution in Figure 03 While in infinite volume one expects the median of the gap to 



scale with the mass p, = Z^ttipcac [28||, in the e— regime p, governed by the IR cutoff of 
the volume, remains finite while mpcAC — > 0. This is clearly the case in our simulations. 
Therefore one does not need a chiral action to study the epsilon regime, though the explicit 
symmetry breaking effects should be small compared to the inverse lattice size. As long as the 
volume is large enough that the NLO relations describe the two-point functions, continuum 
XPT results can be used to analyze Wilson fermion data. Since separating topological 
sectors with Wilson fermions is not always possible, we analyze our data averaged over the 
topological charge. 

For completeness we give the relevant formulas for two degenerate flavors, averaged over 
the topological charge. The isotriplet pseudoscalar meson correlator up to 0(e A ) is 



Lt , , , t 



P 1 F 2 L 3^^ Lt 



bM-r) + O^)) , (3) 



where P(x,t) = ^|A l 75^ is the pseudoscalar density operator and 



Op = |/i(u). (4) 
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with 

p = 1 



2(FL) 2 

the shape factor (Pi = 0.14046 for our symmetric geometry ) and 

u — 2 mEV^ p . 

I\ can be expressed in terms of Bessel functions, h(u) = 8Y'(u)/(uY(u)). It decreases 
smoothly from 2 at u = to 0.68 at u — 10, the largest value we encounter. The function 

Hr) = \\(r-\f-^] (5) 

describes the quadratic time dependence. The pseudoscalar correlator is dominated by E, 
the dependence on F is only through the 0(e 2 ) term. 

The flavor triplet axialvector current correlator at NLO is 

r A (t) = p J d 3 x{A (x)A (0)) 

= £( a » + ?fef A (i;) + < £4 >)- < 6 » 

with Ao(x, t) = ^|A l 7o75'0 and 

K = i« 2 /i(«) ■ (7) 

In our L s = L t case koo = fti/2. The axialvector correlator is dominated by F, the depen- 
dence on E enters only through the combination mTjV. 

The e-expansion formulas are systematic expansions in the parameter 1/(FL) 2 = C(e 2 ), 
but depend on the 0(1) quantity mEV. In our simulation we explore the range rriEV ~ 
0.7 — 5.0. Large values introduce large NLO and NNLO corrections to the correlators, and 
at some point one transitions into the large volume p-regime. Only by examining the fit 
results will we be able to decide what range of mHV values are acceptable in the e-regime. 

The lattice correlators have to be multiplied by the renormalization factors Zp and Z\ 
to obtain the continuum ones in Equations [3] and [H while in the product mY>V the quark 
mass can be expressed in terms of the PCAC mass as m = tjipcacZa/Zp. In our fit we use 
the combinations 

T A = ZlT^ (8) 

and 

m 2 T P = Z>4 CAc rJ att) (9) 
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Figure 5: The pseudoscalar ( red diamonds) and axialvector (blue bursts) lattice correlators and 
the combined fit results at k = 0.1278 and 0.1279 on the 16 data set. The axialvector correlators 
are multiplied by the factor 50 to better match the scale of the pseudoscalar. 

that depend only on F and mEK, and do a combined fit to Equations [8] and M 

The results of the combined fits on the 16 4 lattices are shown in Figures [5] and El where we 
plot both the pseudoscalar and axialvector correlators (the latter is rescaled by a factor 50 to 
match the scale). We use the time slices [5, 11] in the fit. The data are well described by the 
NLO formulas at all four mass values. The results are summarized in Table HT1 where we list 
not only the predicted low-energy parameters but the combination mSV and an estimate 
for m w L as well. We estimate the infinite volume pion mass using the GMOR relation 
m 2 = NfmH/F 2 + 0{m 2 ). In the e-regime one requires m n L <C 1 , though according to the 
analysis of Ref. (4J1, the e- and p-regimes connect smoothly around m n L ~ 2, so values up to 
that level are also acceptable. While at k = 0.1278 both mTV and m^L are somewhat large, 
the fit indicates that the data are described well by the e-regime forms at all k values. The 
predicted low-energy constants, especially F, show a slight drift as m decreases, indicating 
that higher order effects are nevertheless not negligible. Considering that the expansion 
parameter 1/(FL) 2 ps 1.45 is not at all small, this is quite possible. The 0(e 2 ) corrections 
to the pseudoscalar correlator at t = N t /2 in Equation [4] are 34% at k = 0.1278, decreasing 
to 16% at ac = 0.1281, while the constant term a a of the axial correlator in Equation [7] has 



1 1 1 1 


1 1 1 1 


III 


-_ <$> 
- 


tc=0A278 


1 1) 

r 

> J ^'"A —_ 


1 1 1 1 


?^ ~ 

1 1 1 1 


iii 




11 



< 



0.03 

0.02 

0.01 

0.00 
0.03 

0.02 

0.01 

0.00 



: 1 4 1 1 


i i i i 

K—U. 1 2o0 


1 1 1 i 















— -$> — $ — <$> — * 


;E A 




i i i i 


i i i i 


l l l l 









10 



15 



1 1 1 1 

- <> 


1 1 1 1 

_ /c=0.1281 


iii 

<> ~ 


1 1 1 1 


m x m x z 
i i i i 


iii 







10 



15 



Figure 6: Same as Figure E but at K = 0.1280 and 0.1281. 

30-25% corrections in the same range. A volume of (1.85fm) 4 is not large enough to suppress 
finite volume effects. 

Our second data set is 24 4 , (2.77fm) 4 , considerably larger. The 0(e 2 ) corrections are 
reduced to ~ 15% for the axial correlator, though the corrections are still large, 10-25% 
for the pseudoscalar correlator at our mass values. Figures [7] and M show the result of the 
combined fit. Again, we find good agreement for all correlators in the range of 5 < t < 19, 
though we use only time slices [8, 16] in the fit. The statistical errors are under control 
everywhere, though they increase as the reweighting range increases . 

The data points for t < 5 and t > L t /a — 5 do not follow the e-regime ;\TT predictions. 
A natural explanation is that the heavy excitations that couple to the operators die out 
only at t > 5 and influence the correlators at small distances. If that is indeed the case, the 
correlators should show similar behavior on the 16 4 and 24 4 sets. Indeed, at n = 0.12805 and 
Kj = 0.1281, where we have results on both volumes, both the pseudoscalar and axialvector 
correlators are identical within errors for t < 5, showing the same transient behavior. It is 
somewhat puzzling why a recent result using overlap fermions at similar lattice spacing and 
even smaller quark masses see transient behavior in the pseudoscalar channel up to t w 12 
JJJ. It might be due to the small spatial extent (L s = 16) or the asymmetric geometry 
{L t = 32) used in Ref. [IJ], or that the overlap operator is more extended and excited states 
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Figure 7: Same as Figure but on the 24 4 data set at ft = 0.12805 and 0.1281. 



K 


L 




m n L 


F[MeV] 


E^fMeV] 


0.1278 


16 


3.1(2) 


3.14 


90(3) 


256(6) 


0.1279 


16 


2.1(1) 


2.58 


86(4) 


254(6) 


0.1280 


16 


1.4(1) 


2.11 


83(6) 


252(7) 


0.12805 


16 


1.0(1) 


1.78 


82(7) 


250(7) 


0.1281 


16 


0.68(5) 


1.47 


76(10) 


251(7) 


0.12805 


24 


5.2(3) 


2.71 


90(3) 


248(6) 


0.12810 


24 


3.4(2) 


2.19 


89(4) 


250(6) 


0.128125 


24 


2.6(1) 


1.91 


89(6) 


248(6) 


0.12815 


24 


2.3(1) 


1.80 


92(8) 


245(8) 



Table II: Results from the combined fit to the pseudoscalar and axialvector correlators. The values 
of F and S are converted to physical units using ro = 0.49fm. The combination mEF is predicted 
by the fit while for m n L we estimate the infinite volume pion mass from the GMOR relation. 
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Figure 8: Same as Figure [5] but on the 24 4 data set at k = 0.128125 and 0.12815. 
die out slower. 

We also have measurements of the vector current correlator. The data and the fit quality 
are similar to the axialvector we presented above. Since it does not improve the determina- 
tion of the low-energy constants, we do not include it in our analysis. 

One might ask if the data, especially on the 16 4 volumes, could be better fitted with 
NNLO, C(e 6 ) forms. For the pseudoscalar these two-loop results are known in the continuum 



271 ]. and contain two new low-energy constants, L 3 and L 4 . On the lattice additional terms, 



due to the chiral symmetry breaking of the Wilson action, also enter. For the axialvector 
correlator only NLO results are available. Considering the number of unknown parameters, 
it is not obvious that a meaningful fit could be done even if full NNLO formulas were 
available, though it would be very interesting to test. 

We close this section with a combined plot of the low-energy constants obtained on the 
two volumes, as the function of the parameter mLV (Figure [9]). 

As is evident both from Figure [9] and Table HH the different quark mass data on the 
24 4 ensemble are consistent for both low-energy constants and the results for the chiral 
condensate are consistent on the two volumes. F, on the other hand, shows a drift as mSK 
decreases. Without a large volume data point at mEF < 2 we cannot tell if this is due 
to finite volume effects, or signals the breakdown of the e expansion for mTiV > 2. %PT 
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Figure 9: The low-energy constant F and E 1 / 3 as the function of the parameter m'SV, predicted 
by NLO xPT. 



formulas that connect the e and p regimes could help to decide this issue. Since the next- 
to-leading order corrections to F are over 10% on the 16 4 data set, we prefer using the large 
volume data to arrive at our final prediction, 



F = 90(4)MeV, S 1/3 = 248(6)MeV 
Fr = 0.224(10), S 1/3 r = 0.617(15) . 



(10) 



The errors only include the statistical uncertainties. 

Let us finally compare our results to other recent two flavor computations, even though 
direct comparisons are problematic due to different systematic errors. In the p-regime with 
maximally twisted mass fermions, the ETM collaboration gets in the continuum limit Fr 
0.188(2) (7) and S 1 / 3 ^ = 0.597(9) (15) and a compatible number for £ in the e regime 
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29J] 1 . Another e regime computation has been performed by JLQCD with dynamical overlap 



fermions at fixed topology in a L 3 x 2L, L = 1.8fm box [l4|. They get Fr = 0.217(14) and 
S 1 / 3 ^ = 0.596(10). Given the statistical and systematic errors, these results nicely agree 
with our determination. 

Another method of extracting the low-energy constant is by looking at the distribution 
of the lowest eigenvalue of the Dirac operator and comparing it to predictions from random 
matrix theory. To our knowledge, there are two such results with renormalized Nf — 2 



30fl, JLQCD compute r £ 1/3 = 0.624(17) (27) using these methods. Ref. [3l| 



results. In Ref. 

finds r F = 0.213(11) and r £ 1/3 = 0.594(13) using nHYP link dynamical overlap fermions. 
Again, there is good agreement to our findings. 



V. CONCLUSION 

The data presented in this paper has been generated with moderate computer resources. 
This was possible due to the good chiral properties of the action which come at relatively 
low cost due to the simple nHYP smearing procedure, and the effective reweighting that 



1 The number for E is not explicitly given; the quoted errors are obtained using tq = 0.433fm 
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allowed us to lower the quark mass even further. Obviously, there are still shortcomings 
of our analysis. The simulation is done at just one lattice spacing, so we are not able to 
take our results to the continuum limit. However, HYP smearing has improved the scaling 
properties of a variety of actions. We are confident that also here cut-off effects will be small. 

Moreover, setting the scale by r = 0.49fm is not satisfactory. The value of r is not 
known to high accuracy. We could use F as scale parameter. Apart from that, the e regime 
setup makes it problematic to use other widely used scales like the mass of the Q baryon. 

We still are at finite lattice size and e-regime xPT is a slowly converging expansion in 
1/(FL) 2 . Here, our large volume puts us into a good position and the comparison between 
the L/a — 16 and L/a = 24 results shows that the finite volume effects are under control. 
However, statistical errors due to the limited number of gauge configurations are too large 
for a more substantiated claim. 

For the e expansion to be valid, the parameter mSV has to be 0(1). Our data span 
the range 0.7 to 5.2 and might go beyond the validity of the analytical expressions. An 
expansion that connects the e and p regimes would be very useful to control this aspect of 
the calculation. 

Nevertheless our results are encouraging. We find that reweighting works on a fairly large 
volume of (L/a) 4 = 24 4 , L « 2.8fm, and the statistical fluctuations are under control despite 
quark masses as low as 4MeV. Repeating the calculation at a smaller lattice spacing would 
not be prohibitively expensive and could improve on all of the above mentioned issues. 
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